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Abstract 

Let A and B be two finite sets of points witli total cardinality n, the many to 
many point matching with demands and capacities matches each point ai E A 
to at least and at most a'i points in B, and each point bj G S to at least 
{3j and at most /3'j points in A for all 1 < z < s and 1 < i < t- In this paper, 
we present an upper bound for this problem using our new polynomial time 
algorithm. 

Keywords: 

many to many matching, Hungarian method, bipartite graph, points with 
demands and capacities, minimum perfect matching 



1. Introduction 

A matching between two point sets A and B defines a relationship be- 
tween them. The concept of the matching is used in various fields such as 
computational biology Jl[, pattern recognition jij, computer vision jsl, music 
information retrieval jj and computational music theory jsj. A many-to- 
many matching between A and B maps each point of A to at least one point 
of B and vice- versa. Eiter and Mannila j6| reduced the many-to-many match- 
ing problem to the minimum-weight perfect matching problem in a bipartite 
graph and solved it in 0{n^) time. 

Consider two point sets A = {ai, 02, ... , Og} and B = {bi, 62, • • • , bt} with 
s + t = n. Let Da = {ai, tt2, • • • , Q^s} and = /32, • • • , A} denote 
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the demand sets of A and B, respectively. Let Ca = {a'l, a'2, • • • , o's} and 
Cb = (^'2, ■ ■ ■ , be the capacity sets of A and B, respectively. The 
minimum-cost many-to-many matching with demands and capacities, here 
called MMDC matching, is a matching that matches each point Oi & A to 
at least and at most a'i points in B, and each point bj G i? to at least 
I3j and at most 13' j points in A, such that the sum of the matching costs 
is minimized. A function 6{ai,bj) represents the cost of the pairing {ai,bj). 
Note that S{ai, bj) can be zero, positive or negative. The cost of a matching 
is the sum of the costs of the pairings (oj, bj) for all 1 < -i < s and 1 < j <t. 

Schrijver 0] proved that an MMDC matching can be found in strongly 
polynomial time. In this paper, we propose a new 0{n^) algorithm for the 
MMDC matching problem based on the well known approach of Eiter and 
Mannila ^ and provide an upper bound for it. 

2. MATCHING ALGORITHM 

We construct a complete bipartite graph such that by applying the Hun- 
garian method on it the demands and capacity limitations of the elements 
be satisfied. In the following, we explain how our complete bipartite graph 
G is constructed. 

We represent a set of the related nodes using a rectangle, each connection 
between two nodes with a line, and each node with a circle. So a connection 
between two nodes is shown using a line that connects the two corresponding 
circles. A directed line is used to show a branched line such that the input 
line branches into the output lines. A complete connection between two sets 
is a connection where each node of one set is connected to the all nodes of 
the other set. We show a complete connection using a line connecting the 
two corresponding sets. 

Let U T be a bipartition of G, where S = (\J'^^ Ai) U (U^^^ A'i) U 
(Uj=i -^j) ^ (Uj=i ^j) T = Ui=i Bseti. The points of the sets Ai, Bseti, 
and A'i for all 1 < i < s are called the main points, since they are copies of 
the input points. On the other hand, the points of the sets Xj and Wj for 
all 1 < i < t are called the dummy points. All edges (a, b) that their both 
end points are main points, that is a G AiU A'- and b G Bseti for 1 < i < s, 
are called the main edges. 

The Hungarian method computes a perfect matching where each node 
is incident to a unique edge. We aim to find an MMDC matching that in 
which two or more nodes may be mapped to a same node, that is a node 



may be selected more than once. So our constructed graph contains multiple 
copies of each elements to simulate this situation. Let Ai = {an, . . . ,aia^} 
for 1 < z < s be the set of the ai copies of the point Oj. Each set Ai is 
completely connected to the set Bseti = {bu, . . . ,bti} for 1 < z < s. This 
complete connection is shown using a line connecting the two corresponding 
rectangles of Ai and Bseti. We define a function Fj for A^ sets, such that 
Fi{aik,bji) = 5{ai,bj). Each Ai set guarantees that each point ai & A is 
matched to at least elements of B. 

Note that each node of A has a limited capacity, i.e. it must be matched to 
at most a given number of the points of the other set. Each point is copied 
{a'i — ai) times and constitute the A'i set. Let A'^ — {a'n, . . . ,a'i(^a'i-ai)} 
for 1 < i < s. A'i sets guarantee that each point bj E B is matched to at 
least /3j elements of A. Moreover, each set A'i assures that each point ai is 
matched to at most a- elements. Each A'i is completely connected to Bseti, 
where the cost of {a'id, bji) edge is equal to S{ai, bj) for all 1 < < {a'i — ai). 

Assume that all nodes bji for 1 < i < s constitute sets, called Bj. In fact, 
the set Bj is s copies of bj. We use the Wj = {wji, . . . , set to limit 

the number of the points that can be matched to bj E B for 1 < j < t. There 
is a zero wighted complete connection between the nodes of Bj and Wj for 
l<J<t. 

Let Xj — {xji, . . . ,Xj(^^'.^p.)} and 7 = minS{ai,bj) for all 1 < i < 
s and 1 < j < t. Select an arbitrary number 7' such that 7' < 7, there exists 
a 7' weighted complete connection between the nodes of Bj and Xj for all 
I < j <t. Xj sets guarantee that the matching is a minimum cost matching. 

There exists another set that compensates the bipartite graph, called Y. 
The input of the Hungarian algorithm is a complete bipartite graph, i.e. both 
parts of the input bipartite graph have an equal number of points. So we 
should balance the parts of our constructed bipartite graph before using the 
Hungarian algorithm. We have 
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Figure 1: The constructed complete bipartite graph by our algorithm in which ^l^i ct'i > 
, and 

s 

\T\ = \ [jBseti\ = {s*t). 

i=l 

Let |y| = I JZl^iCt'i — X]j=i/9il, we add F to 5 or T depending on the 
values of X)i=i «i and J2]=i Pj- 

Two cases arise: Y^l^i a- < 13 j and Y^l^^ a- > jij- In the first 
case 1 5*1 < \T\, so wc add the Y set to S. There is a zero weighted complete 
connection between Bseti and Y for all 1 < z < s. 

Now assume that Yll=i^'i > Sj=i/^j- this case, the compensator set 
Y is inserted to T. In this case, there is a complete connection between Xj 
and Y that in which the cost of the edges is an arbitrary number 7" with 
y < 7" < 7. Consequently, the first priority of the points of Xj is the points 
of Bj set. Moreover, A'i is completely connected to Y with zero weighted 
edges. The complete bipartite graph that in which Yl!i=i Oi'i > Sj=i f^j 
shown in Figure 1. 

We claim that from a minimum weight perfect matching in G = S'UT, we 
can get an MMDC matching between A and B. Let M be the union of the 
main edges of a minimum weight perfect matching in G. In the following. 



we prove that the weight of M is equal to the cost of an MMDC matching 
between A and B, caUed L. 

Lemma 1. w{M) < c{L). 

Proof. We get from L a perfect matching M' in our complete bipartite 
graph G, such that the cost of the union of the main edges M" be equal the 
cost of L, that is w{M") = c{L). So we want to convert an MMDC matching 
between A and S to a perfect matching in G. In a perfect matching each 
vertex is incident to exactly one edge, so we relate an edge of G to each 
pairing of L as following. 

Let Pi be the number of the points bj G B that are matched to G A in 
the MMDC matching L. It is Obvious that O-i < Pi < a'i. For each pairing 
(ttj, bj) in L, we connect hji to one of unmatched points of Aj, that is a^^ with 
1 < < cuj. Then, depending on the value of Pi two cases arise: 

• Pi = ai- In this situation, we add the zero weighted edges of G that 
connect each a'y G A'i to one of unmatched nodes of Y for all 1 < j < 
{a'i - ai)- 

• Pi > ai- In this case, we need to match the Pi — nodes of A'i with 
the nodes of Bseti. So, for each pairing (a^, bj) of the Pi — ai remaining 
pairings we add an edge of G that connects a'ij to bji. Then, if yet 
there exist nodes of A'i that have not been matched, for each of them 
we select a zero weighted edge of G that connects it to an unmatched 
node of Y, and add it to the matching M'. 

Now, for each Wjk G Wj we add the edge of G that connects it to an 
unmatched node of Bj. The points of Xj are matched to the points of Bj, 
unless no points remain unmatched in Bj. So we first add the edges that 
connect the nodes of Xj to the remaining unmatched nodes of Bj. Then, 
we add the edges that connect the unmatched nodes of Xj, if exists, to the 
unmatched nodes of Y. 

Note that some points of Bseti sets for 1 < i < s may remain un- 
matched. We match the remaining unmatched points of the Bseti sets 
with the unmatched nodes of Y. Recall that this situation arises when 

Since all nodes of G are selected once, M' is a perfect matching. For each 
{tti, bj) G L there is an edge with equal weight in M", so w{M") = c{L). M 



is the union of the main edges incident to the nodes of M" in a minimum 
weight perfect matching in G, that is w{M) < w{M"), so w{M) < c{L). □ 

Lemma 2. w{M) > c{L). 

Proof. From the union of the main edges M of a minimum weight perfect 
matching in G, we get an MMDC matching L' between A and B, such that 
the cost of L' be equal to the cost of M, that is w{M) = c{L'). For each edge 
m e M, a m — {ttik, hji) oi m = {a'ik, bji) then we add the pairing (oj, bj) to 
L'. Otherwise, no pairing is added to L'. 

For each ai E A for 1 < i < s, there exists the set Ai in G with ai 
nodes which is connected only to one set, Bseti. So the nodes of each Ai 
are selected by some nodes of Bseti, that is bji for I < j < t. Hence, each 
ttj e A for 1 < i < s is selected at least times, so the demand of the point 
tti is satisfied. In G there exist a, plus cu^ — ai copies of each point a^, that is 
the points of A plus the points of So the number of the points that are 
matched to each E A is at most a^. 

Consider the sets Bj with I < j <t, recall that Bj = {bji\l < i < s} and 
the points of Wj sets are connected to Bj for 1 < j < t by zero weighted 
edges. Wj is connected only to Bj, so Wj selects s — P'j number of the s 
members of Bj and 13' j points remain unmatched in Bj. Suppose that k 
points of 13' j points in Bj are selected by the points of Ai sets for 1 < i < s, 
so the /3'j — k remaining points of Bj should be selected by the other sets 
that are connected to it. We discuss two cases, depending on the value of k. 

• if k < (3j then {13' j - k) > {(3'j - f3j). Then, Xj selects the P'j - 13 j 
elements of the remaining members of Bj . we have 

{(3'j -k)- {(3'j - (3j) = P'j -k- (3'j + (3j ^pj-k>0, 

so the remaining (3j — k members of Bj are selected by the points of 
A'i sets. Note that k points of bji points for all 1 < i < s are selected 
by the points of Ai sets and (3j — k points of them are selected by A'i 
sets. The demand of the point bj is satisfied, since (3j — k + k — (3j. 

• if k > (3j then {(3'j — k) < {(3'j — (3j) and all the {(3'j — k) remaining 
members of Bj are selected by the points of Xj, so the result is minimum 
cost result. 



The cost of L' is equal to the weight of M, since for each edge of M we 
add a pairing with equal cost to L', so c(L') = c(M). On the other hand, L' 
is a many to many matching that satisfies the demands and capacities of A 
and B. L is an MMDC matching between A and S, so c(L) < c{L'). Thus 
c(L) < c(M). □ 

Theorem 1. Let M be the union of the main edges of a minimum weight 
perfect matching in G, and let L be an MMDC matching between A and B. 
Then, w{M) = c{L). 

Proof. By Lemma 1 and Lemma 2 w{M) < c{L) and w{M) > c{L), respec- 
tively. So w{M) = c{L). □ 
The time-complexity of the Hungarian method is 0{n^), where the num- 
ber of the nodes of the input graph is 0{n) 0|. The number of the nodes of 
our complete bipartite graph is 0{n^), so the complexity of our algorithm is 

3. Conclusion 

We presented an 0{n^) algorithm for the many to many matching with 
demands and capacities. In this version of the many-to-many matching prob- 
lem, the points of the two sets A and B are ra-dimensional points, with n > 1. 
We can limit their dimensions to one or two dimensional. It is expected that 
the complexity of the n-dimensional matching problem will be reduced by ex- 
ploiting the geometric information. So the one and two dimensional version 
of this problem remains open. 
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